Breakdown of the perturbative renormalization group at certain quantum critical 

points 



O 

o 

(N 

1—5 









I 

o 
o 



> 

o 

en 

o 
o 

-I— • 



I 

O 

o 



X 



D. Belitz\ T.R. Kirkpatrick^, and J. RoUbiihleri 

^Department of Physics and Materials Science Institute, University of Oregon, Eugene, OR 97403 

^Institute for Physical Science and Technology, and Department of Physics 

University of Maryland, College Park, MD 20742 

(Dated: February 2, 2008) 

It is shown that the presence of muhiple time scales at a quantum critical point can lead to a 
breakdown of the loop expansion for critical exponents, since coefficients in the expansion diverge. 
Consequently, results obtained from finite-order perturbative renormalization-group treatments may 
be not be an approximation in any sense to the true asymptotic critical behavior. This problem 
manifests itself as a non-renormalizable field theory, or, equivalently, as the presence of a dangerous 
irrelevant variable. The quantum ferromagnetic transition in disordered metals provides an example. 
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The phenomenon known as critical slowing down is 
a dramatic effect near critical points 1]. It occurs be- 
cause the order parameter (OP) relaxation rate r di- 
verges together with the correlation length ^. The dy- 
namical critical exponent z is defined by r ex ^^ . If soft 
modes other than the OP fluctuations couple to the lat- 
ter, one needs to take them into account as well. For in- 
stance, at the liquicWas critical point one needs to con- 
sider shear modes pj, and the quantum ferromagnetic 
transition |^ |3| is influenced by soft particle- hole excita- 
tions i.4|. These additional soft modes are 'generic' in the 
sense that their softness is not related to the critical point 
(e.g., shear modes in a fluid are soft due to momentum 
conservation), and they introduce diverging tinre scales 
in addition to the critical one 5]. In contrast to the 
diverging length scale ^, which is unique in any given di- 
rection, one therefore usually has to deal with multiple 
diverging time scales. In classical statistical mechanics 
this complication remains restricted to the dynamics, as 
the static critical behavior is independent of it, but in 
quantum statistical mechanics, where the statics and the 
dynamics are coupled, the critical behavior of thermody- 
namic quantities is also affected [3- 

The technical description of phase transitions often re- 
lies on the renormalization group (RG). There are two 
different formulations of the RG. The flrst one, Wilson's 
momentum- shell RG [y, |if , integrates out degrees of free- 
dom with successively larger wave numbers in order to 
derive a simpler effective theory valid only at asymptot- 
ically large length scales, i.e., at criticality. The second 
one, often referred to as the field-theoretic method, has 
its roots in high-energy physics. It introduces an artifi- 
cial 'renormalized theory' in order to study the behavior 
of the physical theory under changes of the momentum 
cutoff, i.e., the basic length scale p|. These two formu- 
lations of the RG are believed to be equivalent |9|. An 
important concept in the field-theoretic formulation is 
the 'renormalizability' of the theory, i.e., the possibility 
of constructing a renormalized theory that is finite in the 



limit of an infinite momentum cutoff by absorbing any in- 
finities into a finite number of renormalization constants. 
Renormalizability is crucial in order to ensure scale in- 
dependent RG flow equations and the power-law critical 
behavior that generically follows fronr them. While this 
concept has no obvious analog in the Wilsonian method, 
it is generally believed that lack of renormalizability can 
be remedied by adding additional renormalization con- 
stants, which in the Wilsonian method corresponds to the 
generation, under renormalization, of additional terms 
that were not present in the original action. In either 
method, one usually proceeds by performing a loop ex- 
pansion in powers of a suitable coupling constant. In 
conjunction with an e-expansion about some critical di- 
mensionality, this allows for a perturbatively controlled 
determination of the critical behavior. 

In this Letter we come to the remarkable conclusion 
that the perturbative RG can break down at a quantum 
critical point if there are multiple divergent time scales. 
Specifically, low-order results of an e-expansion for criti- 
cal exponents may be invalid because coefficients of terms 
that are formally of higher order in e diverge. A second 
conclusion is that nonrenormalizability can be related to 
the existence of a dangerous irrelevant variable (DIV) 
[lOJ rather than to the generation of new terms in the 
momentum-shell method. Our results are generic, but 
for the sake of definiteness we illustrate them by using 
the quantum ferromagnetic transition in disordered met- 
als 1^ ^s ^^ example. This transition is known to have 
an upper critical dimension d+ = 4 |lll Il2| . Accordingly, 
there is an effective quartic coupling constant u that is ir- 
relevant with respect to a Gaussian mean-field fixed point 
(FP) in spatial dimensions d > 4, but becomes marginal 
in d = 4. Standard lore 0, Q suggests that a controlled 
e = 4 — d expansion of the critical exponents is possible 
for d < 4. Indeed, we will show below that to one-loop 
order one finds a flow equation 
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u^ + Oiu^), 



(1) 



where £ = \nb with b > the RG length scale factor. 
This suggests a FP value u* = 12e/ll + O(e^), which 
determines the critical exponent 77 via r] = u* /6 + 0{e^), 



agators 



77 = 2e/ll + 0(e2). 



(2a) 



The correlation length exponent v and the dynamical 
exponent z are given in terms of 77, 



z = 4- ?7 , u = 1/(2 - r/). 



(2b) 



These results have nothing in common with the exact 
critial behavior found in Ref. [13, which is not even given 
by power laws. This creates an obvious conundrum, es- 
pecially since the 4 — e expansion is a more standard and 
straightforward method than the one employed in [l3| . 

Equations (|1I2|) . although formally correct, do not con- 
stitute an approximation to the true asymptotic critical 
behavior in any sense. This can first be seen in the flow 
equations at two-loop order, where the loop expansion 
parameter u couples to an irrelevant variable w, 



du 
H 
dw 



1 1 v^ 

_u2 (l-lnw)+0(u4) 

12 32 ^ )^ y /^ 

11 w 



(3a) 



-2w+ ^uw + ^{1- lnw)w + 0{u^). (3b) 



w is irrelevant with a scale dimension —2 -I- 0(e) and a FP 
value w* = 0. It is common that some scaling function 
depends singularly on an irrelevant variable, making the 
latter a DIV with respect to the corresponding observable 
|lO|. Here, however, the coefficients of the loop expansion 
itself depend singularyon w, leading to a breakdown of 
the perturbative RG |l3l |. This phenomenon is rcmini- 
scient of the breakdown of the virial expansion for trans- 
port coefficients in classical statistical mechanics [l^ . In 
either case, a seemingly systematic expansion in powers 
of a coupling constant is ruined by the appearance of 
logarithmic terms. Physically, w reflects the ratio of the 
diffusive time scale of the itinerant electrons, which is 
characterized by a frequency scaling as a wave number 
squared, Vt ~ fc^ and the critical time scale characterized 
byfi-fc^ Hill. 

In the remainder of this Letter we sketch the deriva- 
tion these results from a model field theory which has 
the same structure as the field theory that describes the 
quantum phase transition in disordered itinerant Heisen- 
berg ferromagnets. In addition to demonstrating the ori- 
gin of the breakdown of the loop expansion manifest in 
Eqs. Q, it will show that this breakdown is technically 
related to the field theory being non-renormalizable in 
a way that cannot be cured by introducing additional 
renormalization constants. This non-renormalizability is 
equivalent to the dangerous irrelevancy of w. 

Let us consider a quantum field theory that couples dif- 
fusive excitations g to a (/)^-theory, with Gaussian prop- 
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(4a) 
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The fields q and (j) are functions of a wave vector p and 
a Matsubara frequency ri„ = 2TrTn. The propagators 
in Eqs. (|4a|l and (|4b|l have structures appropriate for 
diffusive particle-hole excitations and paramagnons, re- 
spectively, in disordered metals, t is the distance from 
a critical point for the (/(-excitation, and G, iJ, a, and 
Ci are model parameters. Notice that the propagators 
contain the two time scales mentioned above. This will 
be crucial later. Motivated by the known properties of 
the ferromagnetic problem |ll| we now assume that un- 
der renormalization, t, G, and Gi are not singularly af- 
fected, while a and H acquire momentum and frequency 
dependences, respectively, that are determined by a set 
of coupled integral equations which read, at quantum 
criticality, T = t = 0, and with a momentum cutoff A, 



a{k) = a ~\ — — / duj dp 
3 Jo J\k\ 

3G /"^ 

H{n) = H + ^J dp 



p 



d-1 



|fc| [p2 + H{lu)uj] 



a{p)p'^ + aCifl 



(5a) 



(5b) 



Here G2 is another coupling constant. 

While Eqs. (|4I5|I are sufficient to define our model, it 
is conceptually helpful to consider the structure of a field 
theory that generates these propagators. A Gaussian 
Landau-Ginzburg- Wilson (LGW) action consistent with 
Gaussian propagators given by Eqs. Q has a structure 

Ao = dx (/) {t ~ a dl) (j) + Vr-fi dx 4>q 

+ fdxq{-dl/G + H\n,,\)q. (6) 

Here we consider the fields (f) and q functions of position x 
and Qn , and we use a schematic notation that leaves out 
everything not needed for power counting, including sums 
over Matsubara frequencies and the distinction between 
powers of T and Sl„. The coefficients t, a, G, and H 
are the same as in Eqs. Q), and 71 = i/Gio/G. The 
renormalizations given by Eqs. ((HJ are produced by a 
non-Gaussian term with a structure 



Ant ^VT12 dx (j)q^ 



(7) 



with 72 ~ ^GGija,. It gives rise to renormalizations of 
a and H by means of the diagrams shown in Fig.[T1|,lfij. 

The above model shares all qualitative features with 
the theory of disordered itinerant Heisenbe rg quan tum 
ferromagnets that was considered in Refs. [llj, [l^]. In 





(a) 



FIG. 1: Diagrammatic contributions to a{k) (a) and H{Q) 
(b), respectively. The solid and dashed lines denote q and 
^i-propagators, respectively. 



particular, the integral equations (|3J) are slightly simpli- 
fied versions of those derived and solved in Refs. I12L llTl 
The prefactors in Eqs. (|3J) have been chosen such that the 
one-loop terms quantitatively agree with those obtained 
from the full integral equations. 

Since the perturbation theory is given to all orders by 
the integral equations lO, one can determine the critical 
behavior by solving the latter. This was done in Ref. I12L 
In what follows we instead analyze the theory by stan- 
dard RG techniques [a, Ll l3 order by order in a loop 
expansion. As we will see, this leads to important in- 
sights which one misses by solving the integral equations, 
and which transcend the specific ferromagnetic quantum 
phase transition problem that inspired our model. 

We define the scale dimensions of a length L and fre- 
quency ft or temperature T to be [L] = — 1, [ft] = [T] = 
z, with 2: a so far undetermined dynamical critical expo- 
nent J^] . The scale dimensions of the fields (f> and q we 
define such that a and G are dimensionless. We further 
define dimensionless coupling constants h, ci, and C2 by 



H^A 



2-zr 



Cl ^A 



4-z 



Cl 



C2 = A 



2 + €~Z 



C2, (8) 



and derive RG flow equations for these quantities. 

The properties of our model allow us to determine 
several flow equations and exponent relations exactly. 
Since G is not renormalized, the field q does not carry 
a field renormalization, which means that the diffusive 
propagators retains its Gaussian wave number depen- 
dence. (J3, however, does carry a field renormalization, 
which gives rise to the critical exponent 77 defined by 
a{p ^ 0) ex IpI^^'. The fact that Ci and t are not renor- 
malized implies the exact fiow equations 



dci/de^ (4 



r/)ci , dt/de = (2 - r])t. (9) 



Both z and the correlation length exponent v (defined by 
t oc b^''^) can therefore be expressed in terms of 77 and are 
given by Eq. Ij2b|l . The critical properties of the model 
are thus determined by only one independent exponent, 
e.g., 77. Power counting shows that C2 is not singularly 
renormalized either, so the fiow equation for C2 is also 
known exactly. 



dc2/d£ = -(2-e)c2. 



(10) 



The flow equation for h, as well as 77, need to be con- 
structed in terms of a loop expansion. 

Let us first consider the theory to zero-loop order. 
Since C2 is irrelevant for all d > 2, Eq. H1U|) . there is 



a Gaussian FP with r] = 0, z = i. This is the FP found 
by Hertz Q]. It is not stable, as was first shown in Ref. 
y. In the present context, this can be seen as follows. 

Consider the perturbation theory as generated by iter- 
ating the integral equations Q . For our present purposes 
it suffices to do the integrals in d = 4. With U = C2/aH, 
the perturbation theory then has the structure 

a{k)/a = l + C/F(i)(A/|fc|) + C/2F(2)(A/|fc|) 

+0{U'^), (11a) 
H{n)/H = 1 + UG'^^\A/u;) + U^G^^\A/lu) 

+0{U^), (lib) 

The proper expansion parameter is thus not C2, but 
rather U, or its dimensionless counterpart u — C2/ah, 
which for d < 4 is relevant with respect to the Gaussian 
FP with a scale dimension of e, as can be seen from Eq. 
(jHl). The Gaussian FP is therefore stable only for d > 4, 
despite the irrelevancy of C2, since h acts as a DIV J19IJ . 

To one-loop order, F^^'>{x) = filnx with /i = 1/6, 
and G^^^ {x -^ 00) = qi \nx + 0{x^^) with gi = 3/4. Via 
standard methods 20] this yields Eqs. H1I2|I . 

At two-loop order, F'-^'> {x ^^ 00) — /2,2 In x + 0{lnx) 
with /2,2 = -3/32, and G'-'^^x -* 00) = 32,2 In^ a; + 
.g2,ilnx + 0(1) with 93,2 = -1/16, 32,1 = 1/32. It 
is obvious, in either the momentum-shell or the field- 
theoretic RG methods, that the coefficients of the two- 
loop In -terms must not be independent of the coefficients 
of the one-loop In-terms if the theory is to be renormaliz- 
able (in the latter method), or lead to scale- independent 
flow equations and critical exponents (in either method) . 
A straightforward analysis shows that the condition is 
/2,2 = 32,2 = -figi/2. This condition is obeyed by 52,2, 
but not by /2,2- This problem is germane to the structure 
of the perturbation theory and cannot be cured by intro- 
ducing additional renormalization constants (within the 
field-theoretic method). To see its origin, let us analyze 
the contributions to the In -terms. 

Equation IjSbp for H{fl) in our model has the same 
structure as insertions in a classical theory, since there is 
no frequency integral. The one-loop correction to a{p) 
is linear in ln(A/p) = In A — Inp with p = \p\. Both 
of these terms contribute to the In A contribution to 
H{U), with relative prefactors of 1 and —1/2, respec- 
tively. This structure guarantees that 52,2 = — /i5i/2. 
The one- loop contribution to H{U), on the other hand, is 
linear in ln(A/ri^/''), which leads to contributions to the 
In A term in a{k) with relative prefactors of 1 and —1/4, 
respectively. As a result, one has /2,2 = — (1 — l/4)/igi = 
— 3/151/4. This is a direct consequence of the exis- 
tence of two time scales: If one replaced the ft in the 
integrand in Eq. Ij5b|l by f2^, which would result in a 
model with only one time scale, fi ~ p^, one would get 
/2,2 = 52,2 = -/i.gi/2. 

These observations suggest the following interpreta- 
tion of this failure of the renormalization procedure. Af- 



ter scaling the integration variable u in Eq. H5a() with 
f? jH, the one-loop correction to il{uS) is proportional 
to ln(A4iJ/Ciwp2)^ Writing this as ln(AVp4) _ inw + 

ln(iJp2/Ci), we can write F^"^^ in Eq. (|lla|l as 



F(2)(A/|fc|) = F(2)(A/|fc|)-- 



dp- Inw(p), (12) 
fel P 



where F^^\x -^ oo) = — (lna;)^/16-|-0(lnx), and w{p) ~ 
Hp^/Ci = {h/ci){p/KY. This leads to the two-loop flow 
equations ©J where w — w{p = A) = h/ci. The FP value 
of the irrelevant variable w is w* — 0. Inserted into 
the u-flow equation, this yields an infinite coefficient at 
two- loop order. We also see that the physical origin of 
this problem is the second (diffusive) time scale, which is 
reflected in the scale dimension of w ex /i. 

The following picture now emerges, /i is a DIV since 
the loop expansion is in powers oi u = C2/ah. At one- 
loop order, this makes the Gaussian FP unstable, but 
one has the seemingly sensible FP discussed in connec- 
tion with Eqs. H1I2|I . The two-loop term, however, goes 
as h~^\nh rather than just as 1/h'^, and this leads to 
a breakdown of the loop expansion. If one eliminates 
w by solving Eq. I)3b|l and substituting the result in Eq. 
()8a|l . one obtains scale dependent flow equations and crit- 
ical exponents. As we have seen, this is a consequence 
of the presence of two time scales in conjunction with 
the structure of the perturbation theory for our model. 
It is therefore not possible to analyze this fleld theory 
by means of standard RG methods. At higher order in 
the loop expansion, higher powers of Inw; will appear. 
Clearly, a resummation to all orders is needed in order to 
determine the critical behavior. Indeed, a solution of the 
integral equations, which achieve such a resummation, 
shows that the asymptotic critical behavior consists of 
power laws multiplied by log-normal terms J12l | . 

In summary, we have used the quantum ferromagnetic 
transition in disordered metals to show that one cannot 
necessarily trust low-order perturbative RG results, even 
in the absence of any indications that the perturbative 
FP one has found might not be stable. The mechanism 
we have discussed involves a DIV and is rather general, 
and is thus expected to be operative at other phase tran- 
sitions as well. For instance, our analysis casts doubt on 
a recent one-loop calculation for the clean quantum ferro- 
magnetic problem [21j. We note, however, that the pres- 
ence of multiple time scales, while necessary, is not suffi- 
cient to produce the effect. For instance, at the liquid-gas 
critical point mentioned in the introduction, the analog 
of the irrelevant variable w is not dangerous J23 | . 
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